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Abstract 

Recent results from the PAMELA [l], ATIC [2] , PPB BETS Q and Fermi 
[3] collaborations extend the energy range in the e + ,e~ measurement up to 
unexplored energies in the hundred GeVs range confirming the bump starting 
at about lOGeV already suggested by HEAT [5] and AMS01 data Q- This 
bump can be explained by annihilating dark matter in the context of exotic 
physics, or by nearby astrophysical sources e.g. pulsars. In order to discrim- 
inate between competing models for primary positron production, the study 
of anisotropies complementary to the spectrum determination , shows up as 
new tool to look for the origin of the lepton excess. In this letter we calcu- 
late the contribution to the electron flux given by the collection of all known 
gamma ray pulsars (as listed in the ATNF catalogue) and by annihilating 
dark matter both in case we have a clumpy halo or in case the excess can be 
atributed to a nearby sizeable dark matter clump. We address the problem 
of the electron anisotropy in both scenarios and estimate the prospect that 
a small dipole anisotropy can be found by the Fermi observatory. 
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1. Introduction 

The recent claim from the PAMELA collaboration of an excess in the 
ratio of positrons to electrons plus positrons (positron fraction), seems to 
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support what HEAT and AMS01 data suggested regarding the existence of 
a possible primary cosmic ray electrons source. However some uncertainties 
remain , e.g the possible enhancement of the positron flux due to the rela- 
tively poor knowledge of the propagation parameters set 0]. A large number 
of candidates have been suggested that may be able to reproduce the ob- 
served positron fraction. Among the exotic ones, annihilating dark matter 
may be the most widely studied . Dark Matter is assumed to annihilate in 
the galactic halo and its annihilation products should reach the Earth and be 
measured. As it has been pointed out in previous works [8|, PAMELA data 
seems to accommodate preferably leptonic channels rather than dark matter 
decaying to quark pairs , but typically large boost factors are required to 
adjust the shape to the data. This boost factors are related to an increase of 
the annihilation rate due to e.g. the recently invoked Sommerfeld effect, an- 
nihilation of a non thermal WIMP, or the usual inhomogeneities of the dark 
matter halo ie, clumps, dark matter mini spikes etc... Moreover, to accomo- 
date theory to the existing data, different normalizations are required for the 
different experiments. More conventional scenarios can also reproduce the 
spectral shape. Specifically in this work we analyse the possible contribution 
of nearby astrophysical objects i.e. pulsars , and for leptophilic dark matter. 
In the case of multiple sources, a natural explanation for the spectral shape 
in the ATIC data is found, namely the presence of peaks in the spectrum can 
be reproduced in case we have a convenient selection of pulsars, or in case we 
have a distribution of clumps contributing to the total flux with the required 
luminosity, however it is necessary to fine tune the pulsar parameters to fit 
the peaks present in the ATIC spectrum. On the contrary, the Fermi data 
shows a much flatter spectrum that can be easily reproduced by pulsars with 
standard parametrizations. We calculate the predicted dipole anisotropy 
produced in both scenarios, i.e annihilating dark matter clumps and nearby 
pulsars, and assess the Fermi detection capability of such signature at least 
at the 2 sigma confidence level. 

2. Propagation 

Cosmic ray propagation is a diffusive process due to the random galactic 
magnetic fields. If we denote the number of particles of type % per unit 
volume found in a time t at x with energy E by rii(E, x, t) [GeV~ l cm~^\, the 
evolution equation as initially written by Ginzburg and Syrotvatskii |9( can 
be expressed as: 
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Where V(DVnj) is the diffusion term with T)(E, x, t) [cm 2 / s] the diffusion 
tensor; V(wi») is the convection term with v = v(x) [cm/s] the velocity of the 
galactic wind; -^(brii) represents the energy loss with b(E,x,t) [GeV/s] the 
rate of energy change; — is the fraction of particles lost by radioctive decay 
with a characteristic life-time [s]; f3cN ^2 s<i n s <Ji^ s — ftcNaiUi is the num- 
ber density of particles created or destroyed by spallation processes in the 
propagation due to interactions with the interestellar medium of density N, 
with a characteristic cross section a^ s [cm 2 ] . Finally we have the source term 
Qi(E, x, t) [GeV _1 cm _3 s _1 ] that injects particles of type i into the galaxy. For 
electron and positron propagation, the relevant processes are diffusion, con- 
vection, energy losses and the source term. Propagation can be approached 
in two complementary ways . The first one is to solve the transport equa- 
tion using numerical methods in the same way the package GALPROP 1^ 



does. The second one is to analytically solve the transport equation with a 
set of realistic simplifying assumptions. In this work we use the standard 
GALPROP code to get the positron and electron backgrounds and an ana- 
lytical solution of the propagation equation for the primary positron's flux. 
For primary positron sources, we have used the Green functions formalism as 
described in where the two main processes , diffusion and energy losses 
are considered . The resulting diffusion-loss equation for this process is given 
by: 



d 

—n(E,x,t) = D(E) ■ V 2 n(E, x, t) 



dt 



^(b(E)n{E,x,t)) + Q(E,x,t) (2) 



Where b(E) = — ^| = aE 2 + bE + c ~ aE 2 codifies the energy losses due 
to (a) inverse Compton and synchrotron radiation , (b) bremsstrahlung and 
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(c) ionisation. At the energies we are interested E > 10 GeV the energy loss 
is very well aproximated by the synchrotron losses in the interstellar (ISM) 
magnetic fields and inverse Compton off the CMB and starlight at optical 
and IR frequencies. The energy loss due to these processes is calculated as 



dE 32tt / e 2 \ 2 / B 2 \ ( E 



-c — [w 



dt 9 \mc J \ 8tt J \mc 

= 8xl0- 17 (w + GxlO 11 ^ 1 ^ 2 ^ E 2 « lO-^E^GeVs- 1 ] (3) 

where the energy density w = wqmb + w opt -iR is the energy density of 
the photon background [eV/cm 3 ] with wcmb = 0.25eV/cm 3 and w opt _ IR = 
0.5eV/cm 3 . wb stands for the ISM magnetic field energy density wb = §^ = 
0.6eV/cm 3 for B=5//G. 

For the diffusion coefficient D(E) = D E 5 , three setups MAX, MED and 
MIN can be considered [?} which are consistent with the B/C . 



Table 1: Diffusion Setups 



D [cm 2 / s 


5 


MAX 


1.8 xlO 27 


0.55 


MED 


3.4 xlO 27 


0.70 


MIN 


2.3 xlO 28 


0.46 



We solve the equation for a steady state source (DM positron injection) 
and for a nonstationary source (SNRs and pulsars) assuming free boundary 
conditions. 

3. Astrophysical sources of high energy positrons 

Among the astrophysical objects that populate our Galaxy, many can 
contribute to the positron abundance in cosmic rays, but the required energy 
excludes a large part of them. Following a Hillas argument, the astrophys- 
ical sources able to inject the required order of energy can be found in e.g: 
SNRs and pulsars . Gamma Ray pulsars are expected to produce pairs of 
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electrons and positrons as a result of electromagnetic cascades induced by 
acceleration of electrons in the magnetosphere. The accelerated electrons 
emit curvature radiation i.e. photons, above the pair creation threshold by 
magnetic conversion. Positron injection can be expressed as a power law 
with an exponential cutoff at E c , = E~ a e~ E ^ Ec . Previous works (e.g. 



13(] , [l4( , [l5| , [16(] , [l7| , [18|] , [19] , |20j ) have shown the plausibility of the pul- 
sar scenario as sources of primary cosmic ray electrons. In our analysis, we 
will assume a benchmark model as considered in namely, the "standard 



model" (ST). Although more refined models for electron production in pul 



sars can be considered, some of them (e.g. Harding-Ramaty model jiyj or 
the one devised by Zhang and Cheng [21[ ) produce an e ± output well below 
the observations or a comparable one (e.g. Chi et al. |1 7f] ) . In the latter case, 
the most outstanding pulsars produce similar patterns to those considered 
in the ST model, so we will assume the most simple scenario for positron 
production in pulsars as a benchmark model for the study. 
The ST model assumes that all the rotational energy of the pulsar is lost 
through magnetic dipole radiation. As the rotational energy is given by 
E = IQ 2 /2 (where I~ 10 45 gem 2 is the moment of inertia and Q the spin 
frequency), the spindown power will be E = IQQ. For such a magnetic 
dipole radiator, the energy loss rate can be written as a function of the neu- 
tron star radius R, and a the angle between the dipole axis and rotation 
E = — b 2 R(i n^sm 2 a ^ j e ^ ^ Integrating this expression leads to 

the solution of the rotational velocity of a pulsar where the magnetic dipole 
radiation braking dominates: 

n <«> = (TT^ < 4 > 

where tq = gj^g^ * s a characteristic time taken to be around 10 4 years 
for nominal pulsar parameters. The luminosity of the pulsar can then be 
derived as: 

Integrating this expresion over the pulsar age T, the total energy output 
can be approximated by E out = I j dtQQ w IQ /2. If t/r » 1 then 
Qq ~ resulting in an energy output into electrons of: 
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E out [ST] = f e ±E— 

TO 



(6) 



where the energy budget is determined by the spin down power E, the age 
of the source T and the conversion efficiency into pairs f e ± that is assumed 
to be of a few %. Thus, the e ± source for a single pulsar located at a distance 
r, injecting positrons at time t and energy E can be expressed as: 

Q(E, r,t) = N- L(t)E~ a exp (-E/E c )5(r) (7) 

where N is the normalization factor taken to satisfy the total energy 
release constraint E to t = J dr EQ{E, r)dE and L is the luminosity of 
the source. We have introduced an spectral cutoff at E c =lTeV motivated by 
the ATIC and Hess data. This cutoff will be relevant for young pulsars, as 
old pulsars have a maximal e ± energy below the cutoff due to energy losses. 
In this work we have assumed typical pulsar injection L(r) = L ° n+1 with 

the usual braking index n=3 for magnetic dipole radiation braking, resulting 
in Eq. [3j nonetheless a similar analysis can be conducted for other choices 
i.e. exponential decay luminosity as it is expected from microquasars or e ± 
release from the nebula that surrounds pulsars. 

For this purpose we will consider all the gamma ray pulsars listed in 
the Australian Telescope National Facility (ATNF) [22] pulsar catalogue 0. 
Young pulsars with typical ages lower than 10 4 — 10 5 years are considered to 
be surrounded by the pulsar wind nebula (PWN) or a SNR shell that confines 
the injected electrons before releasing them to the ISM. This has to be taken 
into account when we consider the age of the pulsars that can contribute to 
the electron abundance. In this respect we will consider two collections. In 
the first one we will take a lower bound of 10 4 years constraining our pulsar 
collection to ages between 10 4 and 10 7 years. This constraint allows us to 
accept pulsars like Vela which most probably are still surrounded by the 
PWN. In this scenario, we assume a low conversion efficiency for the young 
pulsars (0(1%)) to take into account the posible confinement of leptons. In 
the second one, we will consider that pulsars with ages lower than 5x1 4 years 
cannot contribute to the bulk of electrons, i.e. we constrain our collection 
to mature pulsars. This introduces an injection delay At between the pulsar 
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birth and injection into the ISM due to the confinement of the electrons in 
the PWN. This delay may be important for young pulsars for wich T ~ At 
but for sufficiently old pulsars we can safely dismiss the delay issue and set 
the injection time at the pulsar age. 

The condition for gamma ray emmision is taken to be as 23] the fraction size 
of the outer gap g = 5.5P 26 / 21 -B 12 4 ^ 7 < 1 in terms of the pulsar period P and 
the pulsar surface magnetic field Byi (in 10 12 G units ) resulting in a collection 
of more than one hundred pulsars, from which 3 lie at a distance < lkpc. 
In the determination of the injection spectral index a, we have to take into 
account the constraints that come from observations of synchrotron radiation 
from SNRs. We could also assume, as in a Harding- Ramaty model 19], that 
the e ± have the same spectral index as gamma rays emmited by pulsars 
(which has been measured by EGRET to be around 1.4-2.2 for energies 0.1 < 
E < lOGeV). Altogether, we shall consider spectral indexes from 1.4 to 2.2 . 
For the sake of simplicity we will assume that our collection of pulsars have all 
the same spectral index a = 1.7 and as an additional simplifying assumption 
we will take a universal set of parameters for the whole collection. In this 
way, neither PAMELA nor ATIC can be reproduced properly, but with a 
little of fine tuning work, the spectral features can be achieved. For example, 
as recently pointed out by |24j , the width of the peak produced by continuous 
injection depends on the characteristic time of the luminosity of the pulsar 
by Atpeak ~ a. Proceeding like this, selecting a number of well known pulsars 



tpeak 



and adjusting their luminosity parameters, the spiky spectral shape of ATIC 
can be achieved without violating the PAMELA constraints. Additionally, 
the energy losses set up the age of the source that produces the peak around 
600 GeV reported by ATIC, provided that it is far (~ lKpc) and bright 
enough (~ 1 order of magnitude in the conversion efficiency). As already 
noted by [l9|, PSR B0355+54 fulfills the requirements. For this pulsar, a 
very large conversion efficiency into pairs is required to account for the ATIC 
peak (0(40%)). On the other hand, the much more statistically significant 
Fermi/LAT data shows a much flatter spectrum. Just taking standard pulsar 
parameters is enough to fit the data without having to resort to very large 
conversion efficiencies so this is the approach we will follow, although it 
should be noted that due to the poor Fermi energy resolution at these energies 
, actually existing peaks could be smoothed down to the observed spectrum, 
proving the issue of normalization of every pulsar to be a subtle one . 
Once we have the positron source we proceed to calculate the number density 
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of positrons by solving the diffusion-loss equation (EqJ2]) for a non stationary 
source . The solution of the equation has been previously derived for a 
burstlike power law injection source with a cutoff Ec ([ill, 3) ■ 



^' r '*) = ^F^3^J (l-^)^«e (8) 
where the distance scale is aproximately, 



DdM ^iD (E)t ^:^- s ( 9, 

as a function of the diffusion index 5 and the maximum energy given by 
the energy losses: E max ~l/at with a ~ 10~ 16 GeV _1 s _1 . 

As the source emits with luminosity L(t), the flux will be given by 
c(>(E,r,t) = L(t')(f)(E,r,t')dt' . We have calculated the pulsar contribution 
to the local electron flux in the case of burstlike injection and continuous 
injection for our collection of young and mature pulsars. Even in the contin- 
uous case , due to the steep of the index of dipolar emission, the injection 
is well aproximated by a burstlike event, being the effect of following an 
approach of continuous injection to soften the posible contributions of the 
individual pulsars. In order to reproduce the spectral features of the Fermi 
data, we assume a MED diffusion scenario with an overall conversion effi- 
ciency of 3%. Due to the relative variability of the collection of pulsars in 
the age/distance parameter space, we adjust the conversion efficiency for a 
few number of objects that show prominent features in the spectrum at the 
considered energies, namely, Geminga, Monogem and J2043+2740 ( Fig. [1] 

)• 

We can consider two sets of pulsars (BM1,BM2) to address the lepton 
confinement in the PWN. Our first set will be made up of mature pulsars 
with ages T> 5xI0 4 years, thus , introducing a delay between the pulsar 
birth and the electron release to the interstellar medium of At ~ 5xl0 4 years. 
In our second set, we consider also younger pulsars, including Vela, with 
ages T> 10 4 years and a delay of At ~ 10 4 years. In this scenario, as 
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young pulsars are still surrounded by the nebula that confines the electrons, 
a smaller efficiency must be called for. The efficiency we will assume for 
this pulsars, including Vela, is 0(1%), and is taken to satisfy the bounds 
imposed by the Hess data. A larger efficiency would imply a non neglible 
contribution above the background of secondary electrons around 2 TeV due 
to the contribution of Vela, which has not been not observed. 

The resulting fluxes are shown ( Figs. [2] and 0) for the considered scenar- 
ios. In order to account for the experimental data in the hundreds of GeV, 
the conversion efficiency into pairs of Monogem and Geminga must be above 
the nominal value of 3% but inside the standard range (l%-30%) considered 



m 



261 ] . This values may change depending on the diffusion setup considered, 
the PWN e ± release delay, pulsar cutoff etc... so it should be considered as 
a single realisation of the multiple possibilities that can reproduce the data 
3- 




Figure 1: Contribution to the electron flux injected by the collection of pulsars considered 
in the BM2 setup. 
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Table 2: Pulsar Parameters 



Pulsar 


d[Kpc] 


T [10 5 years] 


E out BM1/BM2 [10 50 GeV] 


f e ± BM1/BM2 [%) 


Geminga 


0.16 


3.42 


7.4/4.4 


10.0/6.0 


Monogem 


0.29 


1.11 


1.7/2.83 


18.0/30.0 


J2043+2740 


1.13 


12.0 


1.6/1.6 


0.1/0.1 


Vela 


0.29 


0.11 


0.0/0.17 


0.0/1.0 




Figure 2: Mature pulsar contribution to the electron flux for a MED diffusion setup and 
overall efficiency / e ± = 3% . We assume an injection index of a = 1.7 and a cutoff 
E c = ITeV. We consider burstlike injection and the efficiencies have been adjusted to 

f Mono gem = lg%) f Geminga = ^ ^2043+2740 = Q J% ^ reproduce the Fermi and 

PAMELA data, (a): E 3 ■ (e+ + e~) (b): Positron Fraction 



4. Dark Matter source of high energy positrons 

According to the standard ACDM cosmological model, a ~ 22% of the 
energy content of the universe is in the form of cold dark matter (CDM). 
Probably the leading candidates to account for it are weakly interacting 
particles (WIMPS), with the neutralino and the kaluza-klein boson B 1 the 
most extensively studied ones. The relic density of these particles is deter- 
mined by their annihilation cross section. Observations of the CMB and large 
scale structure surveys estimate the relic density constraining the annihila- 
tion cross section up to a canonical value of (av) ~ 3 • 10~ 26 cm 3 s -1 . This 
value sets the rate of production of standard model particles , e.g leptons, 



10 



that can be measured as tracers of dark matter annihilation. 
The signal that results from DM annihilation, depends on the squared DM 
density from the astrophysical side, and on the DM particle mass and cross 
section from the particle physics side. Since the annihilation rate depends 
on the squared density, the presence of dumpiness or substructure imply 
an enhancement of the signal compared to a smooth density distribution. 
The present dark matter structure is considered to have its roots in small 
amplitude quantum fluctuations during inflation. In the accepted "bottom- 
up" hierarchical structure formation, smaller clumps gather together to form 
larger systems , completely determined by the initial power spectrum of the 
primordial fluctuations. Galaxies are thus embedded in large dark matter 
halos that in turn are made up of self-bound substructure or subhalos. The 
mass distribution, abundance and internal structure of clumps is determined 
by means of hi gh-r esolution numerical simulations as the one conducted by 
Diemand et al.|27l. In this work we make use of the mass distribution of 




Figure 3: (a):Mature + young pulsar contribution to the e + + e~ flux for a MED diffusion 
setup , burstlike injection and f e ± = 3% . We assume an injection index of a — 1.7 and 
a cutoff E c = lTet^.The efficiencies have been adjusted to fMonogem = 30%, fctminga = 
6%, /J2043+2740 = 0.1%, fy^ a = 1% to reproduce the Fermi data, (b): Mature pulsar 
contribution to the e + + e~ flux for a MED diffusion setup , continuous injection and 
overall efficiency f e ± — 3% . We assume an injection index of a = 1.7 and a cutoff 
E c = ITeV. The efficiencies have been adjusted to f± e onogem = 15%, fg mmga = 6%, 
/j2043+2740 = 0-1% to reproduce the Fermi data. 
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clumps that results from Diemand's simulation as in Cumberbatch et al. 
281 ] expressed as: 



— a (M/M )- 1 exp[-(M/M cuto// )- 2/3 ] (10) 



d\og(M/M Q ) 

Where the lower mass cutoff is M cuto ff — 8.O3xlO~ 6 M and, inspired by 
Diemand et al. simulation , we take an upper mass cutoff at 10 10 M Q . We 
normalise the distribution such that we have a local clump density (r=r = 
8.5Kpc) of 500 pc~ 3 between 1O _6 M and 1O _5 M . We assume for the spa- 
tially dependent number density a spherically symmetrical Navarro, Frenk& White 
(NFW) profile so the density of clumps is given by 

dn(r,M) p dn D 



d log(M/M ) (r/R) [1 + {r/R)f d log(M/M ) 

Where R=20Kpc and po = 0.86 for a correct normalization. 

For the internal structure within a clump, we adopt a NFW density profile 
that gives a reasonable fit to the lightest clumps in Diemand's simulation and 
reproduces well the outer parts of the large scale subhalos. 

p( r ) = — ; — ; \ r ^° — ; — ; tttt (12) 

^ ' c(r/r 20 o)[l + c(r/r 20 o)] 2 1 1 

where the maximum radius of the clump is defined as the radius r 20 o 
where the density equals 200 times the critical density. We assume a univer- 
sal concentration parameter c for all the clumps, found to lie in Diemands 
simulation within the 1.6 < c < 3.0 range, and a constant density core below 
10- 9 Kpc. 

With all these considerations, the source term can be expressed as: 



l \ rM max 

Q(E,r) = ^P- f(E) f NFW (c,M) 



2 

m DM J Mr, 



dU ^ M) -dlo g (M/M & ) (13) 



dlog(M/M, 



e) 



where f(E) denotes the number of positrons generated per annihilation 
and energy interval (for each branching channel) and f^ FW (c, M) is the in- 
tegrated squared density for each DM clump 
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f NFW (c,M)= I 4vrr'V(0^' (14) 



o 



Once we have the source term we proceed to calculate the flux at the 
Earth by solving the steady state (-^n(E,x) = 0) of EqfJ] as already done 



e.g in [llj] or 12j . The Eqj2] can be greatly simplified translating the energy 
into s(E). 



s (£) = /°°5^dx (15) 
Je 



b(x) 



The equation then reads: 



d 

Aip(s, x) - x) = ~Q(s, x) (16) 

where ip{E,x) = b(E) -n(E,x) and Q(E,x) = ^Q{E,x). 

The solution , is given by the convolution of the Green function over the 



sources: 



(3c [ 

<fi(E,x) = — / G(E,x, s , x )Q(s 0l x )ds dx o (17) 
where the Green function is given by : 



GiExs x ) - b[E{So)] e- l *-* )2/D «" (1 8) 

as a function of the integrating variable sq- Ddiff is the diffusion distance 
scale given by Ddiff = (^(s(E) — s(Eo))) 1 / 2 in terms of the variable s(E) with 
E the observed energy and Eq the energy at injection. 

As previously stated, PAMELA data favours leptophilic DM, ie , can- 
didates which annihilation products are predominantely leptons. The case 



13 



of direct annihilation into electron and positron pairs can provide good fits 
to the ATIC data but it is excluded if we take into account the Fermi re- 
sults due to a sharp drop at the end-point (see e.g. [29]). Therefore we are 
left with annihilations into ^ and r ± which provide a much softer injection 
spectrum and can accomodate both PAMELA and Fermi. For such injection 
spectrums, large annihilation rates are required. The astrophysical boost 
factor, i.e, boost based on the presence of dumpiness, has been proved to be 
insufficient to account for the required normalization to adjust the data. In 
this respect, proposals of particle physics boost factors as velocity dependent 
cross sections, have been recently considered. 

In this work our purpose is to evaluate the expected anisotropy resulting 
from sources that reproduces the observed abundances,so we will assume the 
required normalization to fit the data as a result of any of the considered 
mechanisms that can boost the annihilation rate. Here, we consider a DM 
candidate that annilates into r 1 * 1 with a mass of 3.6 TeV (Fig. Hj). 
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POSITRON FRACTION 

AMS-01 : Aguilar M. et al. 2007 
A HEAT: Barwick S. et al. 1997 

* HEAT: Beatty J.J. et al. 2004 

□ CAPRIGE94: Boezio S. et al. 2000 

• PAMELA 2008 
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(b) 



Figure 4: DM clumpy halo contribution to the electron spectrum as described in the text. 



Annihilation into r + r with a cross section 



4.65xl0" 24 cm 3 /s. The considered DM 



mass Mdm = 3.6 TeV and concentration parameter c=1.6. MED diffusion setup, (a): 
E 3 ■ (e + + e~) (b): Positron Fraction 



We also address the posibility of a nearby DM clump being the responsible 
of the bulk of positrons found in the Fermi, ATIC and PAMELA data. For 
this purpose, we treat here the clump as a point-like object at a given distance 
d contributing to the source term with luminosity L in the same fashion as 
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already done by 3(|. This possibility has been found to be very unlikely 
unless the Sommerfeld effect is at play, but such a source would imprint a 
signature in the electron arrival direction that could constitute a signal of 
dark matter annihilation given the absence of pulsars in the neighbourhood. 
In this sense, the evaluation of such a signal makes the study meaningful. 
The expected flux for a clump of DM annihilating into or ^ with masses 
around 3 TeV and 2 TeV respectively can reproduce the Fermi data with 
the observed drop by Hess , and with slightly different normalizations the 
PAMELA data. As a posibility, we show the flux produced by the annilation 
in the r channel of a DM clump located at 0.9 Kpc with a DM mass of 3.6 
TeV . In this case , as in a pulsar scenario, a MED diffusion setup is apropiate 
, as the MIN and MAX cases show prominent bumps in the spectrum at high 
and low energies respectively, that are not observed in the Fermi data. 




Figure 5: DM clump situated at d=0.9Kpc. MED diffusion setup and luminosity L = 
1.0xl0 9 Mgpc~ 3 . We assume a DM candidate of mass 3.6 TeV that annilates into t + t~ 
with cross section (av) = 3xl0 _25 cm 3 /s (a): E 3 ■ (e + + e - ) (b): Positron Fraction 



5. Anisotropies 

Large scale anisotropies of CR have been measured so far to be less than 



1% (3J|, but it is known that if it were produced by sources with some spatial 



structure, small anisotropies should be present in the arrival directions and 
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could be correlated to the potential sources, whether it be known or not. 
In the case of high energy electrons, very light charged particles, diffusion 
competes with large energy losses resulting in relatively short paths O(Kpc), 
so it is expected that we can use them to sample only sources within a certain 
distance and age. The cosmic ray intensity can, in general, be expanded over 
the celestial sphere in spherical harmonics. At first order, when we have a 
clear directionality, we have a dipole anisotropy as could be the case of a 
single source dominating the spectrum. In this case, the intensity can be 
expressed as 

1(6) = I + 51 cos6 (19) 

where I = l/2(I max + I min ) being this maximum and minimum intensities 
related to a forward-backward measurement. 

The calculation of the anisotropy produced by a single source has been 
carried out by [32| and is given by: 

c 3D | ViV | 

* = --^ (20) 

Where D(E) is de diffusion coefficient and N is the electron number den- 
sity. 

In order to detect a statistically significant anisotropy , at the 2a level, 
5 > 2\/2(N evts )~ l l 2 where N evts is the number of events above a given energy 
threshold ,that is a function of the detector acceptance and the exposition 
time. N evts (E > E th ) = <p(r, E) ■ Acc ■ T exp dE. 

Estimates of the expected anisotro py in the case of a dominant pulsar 



have been previously shown in e.g. [13|.[14J| but they fail to take into account 
the possible effects that could arise from a realistic collection of pulsars, there 
could be e.g. systematic cancellations. We have calculated the anisotropy in 
the case of a collection of pulsars taken from the ATNF catalogue and , for the 
first time, in the framework of dark matter annihilation of clumps thoughout 
the halo, both in the case of a dominating point-like source or a distribution 
of clumps. In the presence of an isotropic background plus a number of 
contributing sources to the total flux, the expected dipole anisotropy is given 
by: 

4>t { } 
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where the index i runs over all the discrete sources that contribute to the 
full dipole. The product (f^r^) represents the projection of the individual 
dipole over the direction of maximum intensity that is energy dependent and 
4>T{E,r,t) is the total flux observed at Earth. The projection can be easily 
calculated taking into account the angular separation 9 of two sources on the 
celestial sphere given by: 

( VA1 + A2 2 \ . . 

6 = arctan . - . — - 22 

\ sin b\ sin 02 + cos Oi cos d\ cos(«2 — ol\) J 

Where Al = cos 2 62 sin 2 (a;2 — cti) , A2 = cos Si sin 62 — sin 5\ cos 62 cos(a2 — 
«i) and the right ascension and declination are denoted by a« and <5;. 
In this work, we evaluate the anisotropy (Fig. [6]) produced by the collection 
of pulsars considered in Fig. [2] (a) for a MED diffusion setup. 
The dipole anisotropy will change direction depending on the energy, but 
the main contributor to the full dipole is shown to be Monogem above a few 
tens of GeV up to the TeV and in second place Geminga. A big contribution 
to the anisotropy (and the positron flux) can be expected also at higher 
energies by younger objects like Vela or CTAl, although it is still unclear if 
such an object can produce a sizeable amount of electrons due to the PWN 
confinement. The full dipole (Fig. [7j) will be given by the projection of the 
individual anisotropies in the direction of the maximum intensity at each 
energy, resulting in a clear signal in the direction of Monogem at energies 
above 20 GeV. The contribution due to Geminga has an addition effect to 
the full anisotropy due to the relative position between it and Monogem, so, 
changes in their individual contributions to the flux , codified for instance 
in the pair conversion efficiency, should not change the anisotropic pattern. 
Contributions from other pulsars don't result in a systematic addition of 
their signal to the full anisotropy into a particular direction, due to their 
spatial distribution, so they constitute a kind of isotropic background. In this 
scenario, measurements of a possible privileged incomimg direction should 
point out an excess in the Monogem/Geminga direction, roughly oposite 
to the direction of the MW. It is possible however that we are observing 
the contribution of some yet undiscovered pulsars that could show up in 
a potential study of anisotropies. In this respect, searches for gamma ray 
sources as the one conducted by Fermi will help to support or disfavour the 
test. 
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In the case of a number of clumps contributing to the full dipole, as in the 
case of Diemands simulation, due to the symmetry of the clump distribution 
we would expect a dipole anisotropy in the direction of the Milky Way (MW) 
center. The expected dipole anisotropy will be given by 

5 = -i / <K£,r)(5 clump (r)n)dr 3 (23) 

<PT J 

where (S c ium P (r)n) = 2 6D . , (r Q — r cosy?) is the projection of the clump 

contribution to the full dipole in the direction of the MW center and r,ip are 
the cylindrical coordinates of the clump distribution. 

If we introduce the clump distribution we have used in the previous sec- 
tion, the expected dipole anisotropy for the direct annihilation channel will 
be : 




Figure 6: Individual dipole anisotropy in the electron + positron spectrum for every pulsar 
considered in Fig. [2 (a) for a MED diffusion setup. Also shown the Fermi sensitivity to 
such an anisotropy at the 2 and 5 a CL in 5 years. 
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£- dl „ g (M/M s )- n ||^/^(c,M) (24) 

It must be noted that in case we have an initial source distribution f(E) 
of electrons as the one that comes from annihilation through the t + t~ chan- 
nel, an equivalent result is obtained after the convolution of the source term 
with the ViV term of the anisotropy definition. The resulting signal (Fig. 
[9]) for the parametrization considered in Fig. H] , would not reach the Fermi 
sensitivity for a 5 years survey as the Fermi spectrum is too flat. On the 
contrary, fits to the much bumpier ATIC data would provide a signature at 




Figure 7: Full dipole anisotropy in the electron + positron spectrum from collection of 
pulsars as considered in Fig. [2] for a MED diffusion setup. Also shown the Fermi sensitivity 
to such an anisotropy at the 2 and 5 a CL in 5 years. 
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the 2a CL. towards the MW center providing a hint to the origin of the 
positron excess. 

We can also consider a DM point source as the responsible of the bulk of 
PAMELA electrons as could be the case of a larg e DM clump 30] or min- 



In this case the expected 



6D(E)d 



be- 



ispikes around an intermediate mass black hole [3_ 
anisotropy is reduced to 5 = ^ £,r ^^" ts °" rce where S pointsource 

ing d the distance to the source (Fig. 
Although the probability of finding such a bright clump in our neighbour- 
hood is rather small, there are some scenarios, with the Sommerfeld effect 
at play, where this probability can be boosted up to a 15%. In this case, 
the anisotropy would exceed the 2a level pointing toward the existence of 
a dominant source. This measurement should also be complemented with 
searches for gamma ray emission to achieve a consistent prediction. 



Dipole anisotropy in a MED diffusion setup 
Fermi (5 years) 95% CL. 
Fermi (5 years) 5a. 




10' 

Energy (GeV) 



Figure 8: Dipole anisotropy in the electron + positron spectrum from the DM point-like 
source as considered in Fig. HlAlso shown the Fermi sensitivity to such an anisotropy at 
the 2 and 5 a CL in 5 years. 
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Figure 9: Dipolc anisotropy in the electron + positron spectrum from a distributed DM 
source considered in Fig. [JJAlso shown the Fermi sensitivity to such an anisotropy at the 
2 and 5 a CL in 5 years. 



6. Conclusions 

A number of posibilities have been proposed as potential candidates to 
account for the positron excess . From the standard astrophysical point of 
view, pulsars seem to be the most promising candidates. Just considering the 
already known gamma-ray pulsars is enough to explain the spectrum for rea- 
sonable model assumptions, but the features in the ATIC data require quite 
a large conversion factor into pairs for the considered model . On the other 
hand, Fermi data shows a much flatter spectrum that can be reproduced 
with a reasonable set of pulsar parameters. If we consider dark matter anni- 
hilation, the contribution from the clumpy halo can reproduce the observed 
patterns once we have the required normalization issue solved by means of 
particle physics boosts on the thermal averaged cross section. One way or 
the other, the measured spectrum can be reproduced and no clear signatures 
can be found to distinguish between a dark matter scenario or pulsars, as the 
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spectral shape seems to lack information of the origin of the electrons. In or- 
der to distinguish between the proposed candidates, the study of anisotropics 
is proposed , as a not performed yet test, that can provide useful information. 
In this work we have analysed the expected anisotropy from a configuration 
of pulsars as in the ATNF catalogue. We have also derived the anisotropy 
in a general dark matter scenario, both in case we have a very bright point 
source and in the case of a clumpy distribution as ilustrated by N-body sim- 
ulations . For the DM point-like source and the pulsar scenarios, the Fermi 
observatory should be able to detect a dipole anisotropy at 2 sigma CL in 
five years at least . On the other hand, we would expect an excess towards 
the MW center in the case of a clumpy halo as the main contributor to the 
primary positron flux, but fits to the Fermi spectrum doesn't provide enough 
events in a 5 years survey. In contrast, the case of the ATIC spectrum would 
open the posibility of direct annihilation into e + e~ implying a much harder 
spectrum. Even annihilations through r ± ,/i ± channels would imprint a sig- 
nature in the electron anisotropy detectable by Fermi that are not expected 
to be observed in view of the the Fermi data. 

Anisotropics have been shown to contribute in a valuable way to disentangle 
the positron excess problem, nonetheless there are still a lot of theoretical 
uncertainties (e.g. the dark matter halo distribution or the mechanism of 
pair production in pulsars). It must be borne in mind that we have not 
taken into account the proper motion of the pulsars or even that of the Dark 
Matter clumps, that could result in an enhancement or suppression of the 
anisotropy. Moreover, the case of a dipole anisotropy toward a pulsar cannot 
exclude a dark matter scenario, as it is possible to have a large Dark Matter 
cloud in the same direction masking the signal. In any case, the precise study 
of electron anisotropies should be conducted together with pulsar surveys to 
help to discriminate between astrophysical and more exotic sources. 
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